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ABSTRACT 

We review a proposed multicomponent model to explain the features of the pulsed emission and spectrum of the Crab Pulsar, on 
the light of the recent detection of pulsed emission above 25 GeV from the MAGIC atmospheric Cerenkov telescope. This model 
explains the evolution of the pulse shape and of the phase-resolved spectra, ranging from the optical/UV to the GeV energy band, on 
the assumption that the observed emission is due to several components, which have spectra modelled as log-parabolic laws. We show 
that the new MAGIC data are well consistent with the prevision of our model. 
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1. Introduction 

Phase resolved energy distributions of pulsars' signals are im- 
portant to obtain information on the radiation processes and ge- 
ometry of the emission regions in the magnetosphere. At y-ray 
energies, at variance with other energy bands, the three bright- 
est sources (Vela, Crab and Geminga) show remarkably similar 
pulse structures with two main peaks at a large phase separa- 
tion ranging from 0.4 to 0.5. This pattern is confirmed by the 
very recent discovery of pulsed emission from PSR J202 1 +365 1 
(Halpern et al. 2008), that has a peak separation of 0.47. The 
double peak structure of the Crab Pulsar (PSR B0531+21) is 
well observed across the entire electromagnetic spectrum, but 
the intensity ratio between the two peaks changes with energy. In 
particular, in the X and soft y-ray ranges the emission of the sec- 
ond peak (P2) becomes higher than the first one (PI), and a rele- 
vant emission from the region between the two peaks (interpeak 
or bridge, Ip) increases like P2. This behaviour is observed up 
to a few MeV, where the pulse shape turns almost sharply to be 
similar, although not equal, to the optical light curve. 

On the basis of a large collection of data, covering the fre- 
quency interval from the optical to the GeV band, we proposed a 
model (Massaro et al. 2006a, hereafter MCCM) able to describe 
the spectral and phase distributions by means of a double two- 
component model. The energy spectra of these components are 
not described by a simple power law, because of the continuous 
spectral steepening towards high energies. We found that a very 
satisfactory model is a parabolic law in a double-logarithmic 
plot, corresponding to a log-normal spectral distribution. 

A useful test to verify the goodness of this multicomponent 
model, and in particular the existence of the two high energy 
components, is the study of the pulse shape at energies higher 
than a few GeV, where P2 is expected to be again the dominant 
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feature as in hard X/soft y rays. Some hints in this direction are 
given by the EGRET pulse profile (Thompson 2004), but statis- 
tics above 5 GeV are so poor that no firm conclusion can be ob- 
tained. Very recently, the MAGIC telescope (e.g. Lorenz 2004) 
has detected pulsed emission from the Crab (Aliu et al. 2008), at 
energies above 25 GeV and with 6.4 standard deviations signif- 
icance, showing the very well established two peak profile, with 
similar amplitudes of the two peaks. 

In this Note we compare the MAGIC data to the predictions 
of the MCCM model, and use them to constrain the high energy 
components. The model is synthetically described in Sect. 2, the 
extension to MAGIC results is presented in Sect. 3 and some 
hypotheses about the origin of these components on the basis 
of some recent works on the high energy emission from young 
pulsars are discussed in Sect. 4. 



2. The multicomponent model of Crab pulsar 

2.1. The two-component model: optical to hard X-rays 

As early presented in Massaro et al. (2000), the spectral and 
phase changes of Crab X-ray pulse shape are well reproduced 
by two components. The first component, called Cq, is assumed 
to have the same pulsed profile observed at optical frequencies, 
while the second component, Cx, is described by an analytical 
model to reproduce the observed pulse profiles. The latter com- 
ponent dominates at the interpeak (Ip) and second peak (P2) 
phase regions. This choice can be justified also from theorical 
models like that proposed by Eastlund et al. (1997), who con- 
sidered the synchrotron emission from marginally clamped elec- 
trons in a shell at the boundary of the closed magnetosphere and 
calculated pulse profile very similar to Co- 
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Fig. 1. Left panel: The ratio between the fluxes of P2 and PI phase regions (PI: -0.06-0.04; P2: 0.32-0.43), compared to the 
predictions of the model. The black data points come from various experiments (Kuiper et al. 2001), while the red point at GeV 
energies has been obtained from MAGIC data. The two extrapolations above 1 GeV correspond to different values of the cut-off 
energy of the Co y spectrum: 15 GeV (solid red line) and 1 1 GeV (dashed blue line). The Cx y cut-off energy is fixed at 15 GeV. Right 
panel: Broadband spectral energy distribution of the total averaged pulse and of the four components: dashed line: Co', dash-dotted 
line: Cx', dot-dot-dashed line: Co y ', dash-dash-dotted line: Cx y - Data points and SED of individual components are the same of the 
MCCM model with the addition of the new MAGIC data (black square, where the shaded area represents the systematic error on 
the spectral fits performed by Aliu et al., 2008). 



By means of spectral fits of BeppoSAX data, we obtained 
that the Co and Cx components were described by a log- 
parabolic spectral law, 

F(E) = KE- {a+hhogE) (1) 

where E is the energy in keV, K is the photon flux at 1 keV 
and the parameter b describes the "curvature" of the log- 
parabola. Best fit estimate of the parameters allowed to derive 
the peak energies of the spectral energy distributions (SED), 
E p = 10 (2 ~ fl)/2 '\ that were found at 12 keV and 178 keV for C 
and Cx, respectively, and b was found the same for both compo- 
nents and equal to 0. 16. 

2.2. The model aty-ray energies 

Observations performed by the COMPTEL and EGRET experi- 
ments onboard the Compton Gamma Ray Observatory (CGRO; 
Kuiper et al. 2001 ; Thompson 2004) provided light curves above 
10 MeV of a good statistical quality, which show that the pulse 
shape returns to be similar to that of Co, i-e. the optical/soft X- 
ray one, although some minor differences are present, for in- 
stance in the shape of P2. At energies higher than ~500 MeV 
there is some indication that the emission from Ip and P2 ap- 
pears to increase, in analogy with the X-ray band. Left panel of 
FigureQ]shows the energy evolution of the P2/P1 flux ratio, com- 
puted in the phase intervals defined by Kuiper et al. (2001) that 
are reported in the caption. 

To model the y-ray emission and the change of the pulse 
shape, MCCM assumed that there are two more, high-energy 
spectral components, Co y and Cx y , with the same pulse shape 
of the corresponding lower-energy components and with spec- 
tral distributions also given by Eq. (1). This extended model has 
six new parameters, i.e. the peak energies, curvatures and nor- 
malizations of the Coy and Cx y components, that should be de- 
termined by data fitting. Given the statistical quality of CGRO 
data, the resulting estimates are much more uncertain than in 



the X-rays. Therefore, we assumed that the curvature parame- 
ters were equal to the Co and Cx ones (b = 0.16), and adjusted 
the normalizations and peak energies to reproduce the observed 
total (phase-averaged) spectrum. Peak energies of Co y and Cx y 
were found to lie around 300 MeV and 2 GeV, respectively. To 
be consistent with the upper limits to the TeV pulsed emission 
(e.g. Lessard et al. 2000) we added also an exponential cutoff to 
both Co y and Cx y , at the energy E c = 15 GeV. This value was not 
compelled by observational or physical reasons, but was simply 
a guess to take into account the upper limits. The MCCM model 
was then able to reproduce both the broadband energy spectrum 
of the total pulse (Figure Q] right panel), and the spectra of the 
PI, Ip and P2 phase intervals. In the left panel of Figure Q] we 
plotted also the fitted and extrapolated P2/P1 ratio for two values 
of the Co y cut-off energy, i.e. 1 1 and 15 GeV. This ratio depends 
on the normalizations and shapes of the Co and Cx components 
in the proper phase ranges, and this is not immediately appar- 
ent in the phase averaged spectrum of Figure Q] (see MCCM for 
details). 



3. The pulse shape and spectrum in the GeV band 

The pulse shape observed above 25 GeV by the MAGIC tele- 
scope (Aliu et al. 2008) is shown in Figure|2] where it is normal- 
ized to the counts of the first peak. Despite the statistical quality 
of data, is well apparent that P2 has an amplitude comparable 
and possibly higher than PI, in agreement with the trend barely 
apparent from the EGRET data (Thompson 2004). 

In Figure |2] it is also plotted the pulse shape obtained from 
our model, with the same parameters' values used in MCCM, for 
an energy of 25 GeV. The cut-off energy (15 GeV) is the same for 
both the Co y and Cx y components. The amplitude of the peaks 
is very well reconstructed; the pulse shape in the MAGIC band 
predicted by our model has small variations with energy, with 
the height of P2 that varies by about 15% between 25 and 60 
GeV. 
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The spectral distribution in the MAGIC range is still poorly 
determined, mainly because of the uncertainties in its energy 
scale. Aliu et al. (2008) combined the very high energy MAGIC 
data with the COMPTEL and EGRET points and evaluated that 
a possibile exponential cut-off can be at E c = 17.7 + 2.8 + 5.0 
GeV, thus compatible to the value assumed in the MCCM model. 
However, the CGRO-MAGIC data are compatible with both an 
exponential or a super-exponential shape, with a preference for 
the former. We reported the flux estimate in our SED model 
(FigureQ] right panel) and found that it is in a substantial agree- 
ment with the extrapolation which was introduced to take prop- 
erly into account the EGRET data. Finally, note that in the 
MCCM model the flux above ~10 GeV is mainly due to the 
Cx y component (Figure 2). This is not in contrast with the ev- 
idence of PI in the pulse profile, because its flux is concentrated 
in a quite narrow phase interval, and Cx y has however a wide 
pedestal extending below PI (Figure |5J. 

The good matching of the predicted pulse shape with the 
MAGIC one implies that the cut-off energy should be similar 
for Co y and Cx y , because a lower cut-off for Coy would imply 
an higher P2/P1 ratio (as shown in Fig. [TJ and therefore a higher 
amplitude for the second peak. 

One of the assumptions of the MCCM model was that the 
phase distributions of the various components do not change 
with the energy. It was essentially motivated by the need to have 
a low number of parameters and supported by the satisfactory 
agreement of the computed pulse profiles with data. However, in 
Massaro et al. (2000) it was already shown that the introduction 
of a mild energy dependence of two parameters of the analytical 
description of Cx would improve the pulse profile modelling, 
particularly in the Ip region. Our extrapolation at 25 GeV pre- 
dicts a higher Ip than that given by the MAGIC data, especially 
in the leading wing of P2. We expect a number of detected events 
in the Ip phase interval, in excess to the off-pulse level, around 
2600, whereas the measured excess is ~900 counts, but still com- 
patible with the off-pulse emission, due to the large uncertainty. 
Note that the S/N ratio does not allow to reach any firm conclu- 
sion: there is, for example, a dip at phase 0.3 (another one is at 
phase -0.15), well below the mean off-pulse level, that reduces 
the Ip content. The assumption that the shape of Cx y is equal 
to that of Cx could not be exactly verified at such different en- 
ergies. They could be different, with a shallower profile for the 
high-energy counterpart, as already suggested in MCCM upon 
considering the Ip spectrum. 

4. Discussion 

On the basis of a large collection of data, from the optical band 
to y-ray energies, we developed a model (MCCM) that describes 
with a good accuracy the spectral and phase distributions. This 
model has essentially an heuristic value, being aimed to establish 
a consistent scenario for the development of a realistic physical 
model. The recent detection of the pulsed emission from Crab at 
energies above 25 GeV by the MAGIC collaboration is a relevant 
result that completes our knowledge on the spectral properties 
of this important source. We have shown here that these results 
are in a good agreement with the expectations from the MCCM 
model and confirm its validity. 

It would be interesting to investigate theoretically the phys- 
ical plausibility for the presence of two couples of emission 
components producing the observed spectra and phase distri- 
butions. A possibility to be further developed arise from the 
similar curvature parameter for the X-ray and y-ray compo- 
nents, that suggests a synchrotron-self-Compton (SSC) mech- 
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Fig. 2. MAGIC observed pulse shape above 25 GeV, normalized 
at the count rate of the first peak, compared with the prediction 
of MCCM (solid line). 

anism for the emission (Morini 1983, Cheng & Wei 1995). A 
stream of relativistic electrons moving along the magnetic field 
lines would undergo a Compton upscattering on the synchrotron 
photon field. A different location and azimuthal distribution in 
the magnetosphere would explain the different pulse shape for 
the Co and Cx components. An alternative hypothesis is that 
the y-ray components are photons emitted by primary electrons 
via the curvature radiation mechanism, and the softer ones are 
the synchrotron emission from secondary pairs. However, in this 
case one would expect a more pronounced curvature of the log- 
parabolic spectra for the latter components. In fact, if the elec- 
tron spectrum is a log-parabola with a curvature parameter r, in 
the ^-approximation the synchrotron curvature b would be equal 
to r/4 (Massaro et al., 2006b), while for curvature radiation this 
parameter would be r/9, because of the dependence on the elec- 
trons' Lorentz factor as y 3 instead of y 2 for the synchrotron pro- 
cess. About the log-parabolic (or log-normal) energy distribution 
of electrons, we recall that it can originate by stochastic acceler- 
ation processes or when the acceleration probability is energy- 
dependent. 

Several models have been proposed to explain the high- 
energy emission from pulsars, both in the outer gap (Cheng et 
al., 1986, 2000) and in the polar/slot gap (Muslimov & Harding 
2003, 2004) frameworks. A common feature of these models is 
that the broadband emission comes from components originat- 
ing from different physical processes. 

Takata & Chang (2007) developed a 3D outer gap model 
based on the 2D analytical solution of the accelerating field 
and particle motion by Takata et al. (2004, 2006) and Hirotani 
(2006). In this model the X-ray emission of PI is due to two 
separate components, with curved, roughly log-parabolic, spec- 
tra and originating by synchrotron emission of secondary pairs 
in different regions of the outer gap, i.e. below and beyond the 
null-charge surface. Their explanation of the y-ray spectrum in- 
vokes other two components due to Inverse Compton scattering 
of secondary pairs and curvature radiation from primary elec- 
trons. Their modelling of P2 and Ip, unlike the MCCM model, 
has different weights for these two components. The increase of 
P2 and Ip with respect to PI in this model is due to the fact that 
the emission comes from regions in the lower magnetosphere, 
where the high magnetic field produces harder spectra. 

More recently, Harding et al. (2008) performed a 3D simula- 
tion of the emission from the Crab pulsar, assuming a Slot Gap 
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accelerator, in which the emission is due to two distribution of 
particles, primary electrons and electron-positron pairs. The X- 
ray emission is due to the synchrotron radiation from the pairs, 
while the y-rays are curvature radiation and the synchrotron ra- 
diation resulting from the resonant cyclotron absorption of radio 
beam photons by the primary particles. 

Future observations, in particular with the Fermi Gamma-ray 
Space Telescope, will be very useful to confirm the spectral cut- 
off observed by MAGIC, to improve the estimate of E c , and to 
bridge the gap between the 100-500 MeV EGRET lightcurves 
and the >25 GeV MAGIC one. It will be possible to follow in 
detail the evolution of the P2/P1 flux ratio and to verify with a 
much higher statistics whether the Ip emission is actually related 
to that of P2. 

Another future interesting test will come from phase- 
resolved X-ray polarimetry, where MCCM expect that the hard 
X-rays polarization properties of P2 should become increasingly 
similar to those of Ip, because of the higher contribution of the 
Cx component. 
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